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INTRODUCTION 


A  test  Droblem  was  reauired  for  the  evaluation  of  the  computer  aided 
design  (CAD)  system,  which  is  currently  being  purchased  by  Benet  Weaoons 
Laboratory.  We  needed  a  relatively  simple. stress  analysis  problem  which 
could  easily  be  done  on  a  CAO  system,  by  the  potential  contractor,  to 
demonstrate  a  particular  system.  The  geometry  selected  was  a  simple  thin 
cylinder  of  mean  radius  1.0  subjected  to  an  axial  tension  load.  The  load  is 
distributed  at  the  ends  of  the  cylinder  to  produce  a  uniform  unit  (1.0)  nomi 
nal  stress.  The  cylinder  is  3  PI  (9.4278)  long  and  contains  a  single  per¬ 
foration  at  mid-length,  in  the  form  of  a  circular  hole  of  radius  0.15.  The 
problem  is  to  find  the  stress  concentration  factor  (S.C.F.)  for  the  hole  whe 
the  thickness  of  the  cylinder  is  0.0143525.  The  most  desirable  solution  is 
for  a  shell  formulation  with  membrane  and  bending  stiffness  terms. 

Along  with  the  S.C.F. ,  it  is  also  desirable  to  generate  the  following 
graphic  output! 

1.  A  deformed  structure  plot. 

2.  A  maximum  principal  stress  contour  plot  on  the  inside  surface. 

3.  A  maximum  principal  stress  contour  plot  on  the  outside  surface. 

If  necessary,  the  problem  may  be  reduced  to  a  membrane  problem  with  only  one 


contour  plot. 


REFERENCE  SOLUTIONS 


There  are  two  readily  available  solutions  to  the  problem  in  standard 
reference  books  (refs  1,2). 

The  solution  in  reference  1  is  a  shell  solution  in  which  the  S.C.F.  is 
plotted  as  a  function  of  a  shell  stiffness  parameter  beta.  In  contrast,  the 
solution  in  reference  2  is  for  a  thick  cylinder  in  which  the  S.C.F.  is  given 
as  an  equation  using  the  internal  and  external  diameters  along  with  the  radius 
of  the  •'ole.  The  rather  strange  thickness  of  0.0145325  comes  from  the  desire 
to  read  the  curves  given  in  Reference  1  at  beta  =  0.80,  which  results  in  an 
S.C.F.  of  4.60.  If  it  is  necessary  to  use  a  membrane  solution  without  bending 
terms,  the  S.C.F.  is  reduced  to  4.00. 

FINITE  ELEMENT  SOLUTION 

The  full  shell  solution  has  been  done  on  Abaqus  at  three  different 
thicknesses  (Figures  1,  2,  5,  and  8).  This  provided  a  reference  solution  to 
aid  the  evalution  committee.  The  S.C.F.  results  are  shown  in  Table  I  below, 
and  Figures  3,  4,  6,  7,  9,  and  10  in  this  report. 


TABLE  I.  STRESS  CONCENTRATION  FACTORS  FROM  SEVERAL  SOURCES 


Thickness 

S.C.F. 
Reference  1 

S.C.F. 
Reference  2 

S.C.F. 

Abaqus 

Percent 

Error 

A 

0.00143525 

8.50 

3.183* 

7 . 20 

-15.3 

B 

0.0943525 

4.60 

3 . 187* 

4.31 

-  6.3 

C 

0.143525 

3 . 25* 

3 . 227 

3.06 

-  5.2  ! 

*These  solutions  are  outside  of  the  published  limits  in  the  reference. 


I  ’lR.  E.  Peterson,  Stress  Concentration  Factors,  John  Wiley  and  Sons,  New  York, 

1974,  p.  153.  . . 

2R.  J.  Roark  and  C.  W.  Young,  Formulas  for  Stress  and  Strain,  McGraw-Hill, 
Fifth  Edition,  1965,  p.  601 . 

t 
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DISCUSSION  OF  THE  PROBLEM 

There  are  two  small  problems  that  must  be  pointed  out  as  they  could  have 
an  impact  on  the  evaluation  of  other  solutions.  First,  the  grid  used  for  this 
analysis  did  not  concentrate  enough  degrees  of  freedom  near  the  hole.  This 
results  in  a  very  poor  deformed  shape  for  solution  'A1  with  its  rather  large 
error.  In  other  words,  the  ‘near  field'  region  of  the  hole  was  larger  than 
was  assumed,  and  element  size  was  increased  too  fast.  In  solutions  'B1  and 
'C',  the  stiffer  shell  has  reduced  the  magnitude  of  this  effect.  The  second 
problem  encountered  was  that  Abaqus  did  not  calculate  the  surface  normal 
direction  accurately  for  all  nodal  points  and  this  data  had  to  be  input 
independently.  This  was  a  rather  large  error  and  until  it  was  corrected, 
solution  ’ B '  was  giving  an  S.C.F.  of  20.2.  Both  of  these  effects  could  be 
corrected  by  the  rather  simple  solution  of  using  a  large  number  of  elements. 
This  method  becomes  very  easy  on  a  computer  graphics  system.  While  a  general 
increase  in  element  density  works  well  on  small  problems,  it  is  not  always 
possible  because  of  the  increase  in  solution  time. 

The  plots  are  taken  directly  from  the  Abaqus  graphics  package  and  are 
plots  of  the  undeformed  grid,  the  deformed  grid,  or  maximum  principal  stress 
contours.  The  grid  is  located  at  the  mean  radius  of  the  cyinder,  but  the 
stress  contours  are  on  either  side  of  the  cylinder  wall.  Contour  plots  for 
section  point  one  are  on  the  inside  of  the  cylinder  and  for  section  point 
three,  they  are  on  the  outside. 
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CHELL  SOLUTION  (A) 


STEP  1  INCREMENT  1 


ABAQUS  VERS  I  ON 


Figure  3.  Near  field  maximum  principal  stress  plot  for  the  inner 
surface  of  solution  A. 
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SHELL  SOLUTION  (A) 


STEP  1  INCREMENT  t 


AflAQUS  version 


Figure  4.  Near  field  maximum  principal  stress  plot  for  the  outer 
surface  of  solution  A. 
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SHELL  SOLUTION  (B) 
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Figure  5.  Deformed  grid  for  solution  8. 
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SHELL  SOLUTION  (B) 


STEP  1  INCREMENT  1 


ABACUS  VERSION  4-5-154- 


Figure  6. 


Far  field  maximum  principal  stress  plot  for  the  inner 
surface  of  solution  B. 
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SHELL  SOLUTION  (B) 
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Figure  7.  Near  field  maximum  principal  stress  plot  for  the  inner 
surface  of  solution  B. 
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SHELL  SOLUTION  (C) 

STEP  1  INCREMENT  t  AflAQUS  VERSION 


Figure  8.  Deformed  grid  for  solution  C. 
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SHELL  SOLUTION  (C) 


STEP  1  INCREMENT  1 


ABACUS  VERSION  4-5-15L 


Figure  9.  Near  field  maximum  principal  stress  plot  for  the  inner 
surface  of  solution  C. 
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SHELL  SOLUTION  (C) 


STEP  1  I  NCR  EVENT  1 


ABACUS  VERSION  4-5-154 


Figure  10.  Near  field  maximum  principal  stress  plot  for  the  outer 
surface  of  solution  C. 
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